Ground state phases of the Half-Filled One-Dimensional Extended Hubbard Model 
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Using quantum Monte Carlo simulations, results of a strong-coupling expansion, and Luttinger 
liquid theory, we determine quantitatively the ground state phase diagram of the one-dimensional 
extended Hubbard model with on-site and nearest-neighbor repulsions U and V. We show that spin 
frustration stabilizes a bond-ordered (dimerized) state for U » V/2 up to U/t ~ 9, where t is the 
nearest-neighbor hopping. The transition from the dimerized state to the staggered charge-density- 
wave state for large V/U is continuous for U J$ 5.5 and first-order for higher U. 
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The one-dimensional Hubbard model, which describes 
electrons on a tight-binding chain with single-particle 
hopping matrix element t and on-site repulsion U, has 
a charge-excitation gap for any U > at half- rilling 
In the spin sector, the low-energy spectrum maps onto 
that of the S =1/2 Heisenberg chain; the spin coupling 
J = 4t 2 /U for U — > oo. The spin spectrum is therefore 
gapless and the spin-spin correlations decay with distance 
r as (— l) r /r 0- Hence, the ground state is a quantum 
critical staggered spin-density- wave (SDW). In the sim- 
plest extended Hubbard model, a nearest-neighbor repul- 
sion V is also included. The Hamiltonian is, in standard 
notation and with t = 1 hereafter, 

H = -* X ^2( c U+i c ^ + 4,i c <M+i) 

<7=T4 i 
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The low-energy properties for V < U/2 are similar 
to those at V = 0. For higher V the ground state 
is a staggered charge-density- wave (CDW), where both 
the charge and spin excitations are gapped. The tran- 
sition between the critical SDW and the long-range- 
ordered CDW has been the subject of numerous studies 
|!SI!ll0SHHmHIIiO3. Until recently, it 
was believed that the SDW-CDW transition occurs for 
all U > at V > U/2 and that it is continuous for small 
U 5) and first-order for larger U. However, based 
on a study of excitation spectra of small chains, Naka- 
mura argued that there is also a bond-order- wave (BOW) 
phase [lfj , where the ground state has a staggered mod- 
ulation of the kinetic energy density (dimerization), in a 
narrow region between the SDW and CDW phases for U 
smaller than the value at which the transition changes 
to first order. Previous studies 0, S H had indi- 
cated an SDW state in this region. Nakamura's BOW- 
CDW boundary coincides with the previously determined 
SDW-CDW boundary. The presence of dimerization and 



the accompanying spin gap were subsequently confirmed 
using quantum Monte Carlo (QMC) simulations |TlLll2^. 
The BOW phase now also has a weak-coupling theory 

The existence of an extended BOW phase has re- 
cently been disputed. Jeckelmann argued, on the basis 
of density-matrix-renormalization-group (DMRG) calcu- 
lations, that the BOW exists only on a short se gme nt of 
the first-order part of the SDW-CDW boundary [jj, i.e., 
that the transition always is SDW-CDW and that BOW 
order is only induced on part of the coexistence curve. 
However, QMC calculations demonstrate the existence 
of BOW order well away from the phase boundary [l2^ . 

Although several studie s ag ree on the existence of an 
extended BOW phase [H El H G3 , the shape of this 
phase in the (U, V) plane has not yet been reliably de- 
termined. The system sizes used in the exact diago- 
nalization study jlfj were too small for converging the 
SDW-BOW boundary (i.e., the spin gap tr ansition) for 
U > 4. In the previous QMC studies [Illll2j . the empha- 
sis was on verifying the presence of BOW order and the 
phase transitions for U ~ 4. In this Letter, we present 
the complete phase diagram. Taking advantage of recent 
QMC algorithm developments — stochastic series expan- 
sion with directed-loop updates |l5| in combination with 
the quantum generalization [rl| of the parallel tempering 
method — we have carried out high-precision, large- 
chain (up to L = 1024) calculations for sufficiently high 
U (< 12) to locate the point at which the BOW order 
vanishes. In agreement with Ref. we nn d that the 
BOW exists also above the U at which the transition to 
the CDW state becomes first-order. However, long-range 
BOW order exists also below this point, and hence the 
point at which the nature of the transition changes from 
continuous to first- order is on the BOW-CDW boundary. 

The phase diagram we find here is qualitatively sim- 
ilar to that obtained in a 4th order strong-coupling ex- 
pansion, where the transition to the CDW state is de- 
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FIG. 1: QMC and strong-coupling phase diagram. The BOW 
is located between the SDW-BOW and BOW-CDW curves. 

termined by comparing the energies of the large- U CDW 
state and the effective spin model including nearest- and 
next-nearest-neighbor interactions J and J' 8]. The 
BOW phase corresponds to the spontaneously dimerized 
phase of the spin chain, i.e., J' /J > 0.241 17]. In Fig. 1, 
we compare our QMC phase boundaries with the strong- 
coupling result; the procedures giving the QMC bound- 
aries will be discussed below. We will show that the sys- 
tem is a Luther-Emery liquid on the continuous BOW- 
CDW boundary, i.e., the charge gap vanishes and the 
spin gap remains open. Evidence supporting this type 
of transition was also presented in Ref. 11]. Here we 
will further argue that the change to a first-order tran- 
sition corresponds to the Luttinger charge exponent K p 
reaching the value 1/4. 

We extract the SDW-BOW and BOW-CDW bound- 
aries using the charge and spin exponents K p and K a . 
If there is a spin or charge gap, the corresponding ex- 
ponent vanishes. Otherwise the equal-time correlation 
function C p (r) ~ r -(*«x C a {r) ~ r -^ 1+K "\ If 
non-zero, the spin exponent K a = 1 as a consequence 
of spin-rotation invariance |18|. On periodic chains the 
exponents can be conveniently extracted from the static 
structure factors S Pt(7 (q) [lflj . 

SpAQ) ^E^'fK'i^lKi^)). (2) 
j,k 

in the limit q —> + : 

K p>a = S p>a (qi)/qi, qi=2ir/L, L -» oo. (3) 

If there indeed are three successive phases, SDW-BOW- 
CDW, as V is increased at a fixed value of U, then the 
spin exponent K a = 1 in the SDW phase and K a = 
everywhere else. The charge exponent K p = every- 
where, except exactly at the BOW-CDW transition point 
if this is a continuous quantum phase transition (i.e., if 
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FIG. 2: Long-wavelength charge (left panels) and spin (right 
panels) structure factors vs V for U — 4 (top), 6 (middle), 
and 8 (bottom). The system sizes are indicated in the low- 
right panel. The U — 4 inset shows the dependence on the 
inverse lattice size at V = 2.10. 

the charge gap vanishes). In contrast, if the transition 
is first-order, then K p = also on the phase boundary. 
Using the relation J3J for a finite system, any discontinu- 
ities will naturally be smoothed, and one can only expect 
to observe Sp <a (qi)/qi developing sharp features as L is 
increased. In Fig. [5] we show results demonstrating this 
for several different system sizes at U = 4, 6, and 8. 

Looking first at the charge exponent, if K p > on the 
BOW-CDW boundary and K p = elsewhere, then one 
can expect a peak developing in S p {q\)/qi versus V. The 
peak position corresponds to the critical V, and the peak 
height should converge to K p . If the transition is first- 
order, S p (qi)/qi should converge to zero for all V, but one 
can still expect some structure at the phase boundary 
for finite L as the nature of the ground state changes. 
In Fig. |3 for U = 4 and 6 the development of sharp 
peaks is apparent. For [7 = 4 the peak-height converges 
to a non-zero value, implying a continuous transition at 
V w 2.160 with K p w 0.43. For U = 6 the convergence to 
a value > is not clear, but the transition point is given 
accurately by the peak-position, which shows very little 
size-dependence. It has been shown previously that the 
transition is first-order for U = 6 0, [n| , and the peak 
should therefore in fact converge to zero. The rather slow 
decay reflects the proximity to the point at which the 
transition becomes continuous. For U = 8 the peak does 
not sharpen, but instead a step develops at the critical 
V. The whole curve converges to as L — > oo. The 
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transition is hence strongly first-order in this case, in 
agreement with previous calculations. As seen in Fig. ^ 
and as observed already by Hirsch |(|, the locations of 
the U = 6 and 8 critical points agree very well with the 
strong-coupling expansion [2(J . 

In the SDW phase, one cannot expect to easily find 
S a (qi)/qi — » 1 exactly, due to logarithmic corrections 
that affect various quantities strongly even for very long 
chains [2ll E3|. However, the log-corrections are known 
to vanish in the frustrated J — J' spin chain at its 
dimerization transition |23|, and hence, since the SDW- 
BOW transition should be of the same nature, the log- 
corrections should vanish here as well. The transition at 
fixed U should therefore be signaled by S a {q\)/q\ cross- 
ing 1 from above as V is increased. Because of the van- 
ishing log-corrections at the transition, the crossing point 
with 1 should not move significantly as L is increased, but 
the drop below 1 should become increasingly sharp, and 
eventually S a (qi)/qi should approach inside the BOW 
phase. This method was used in Ref. E3 and gave a 
slightly higher critical V for the SDW-BOW transition at 
U = 4 than the exact diagonalization ^} . We now have 
results for a wider range of couplings. The results shown 
in Fig. |2 are in accord with the above discussion for all 
three [/-values; S a {qi)/q\ crosses 1 at a T^-point which 
does not move visibly between L = 64 and L = 256. For 
larger V, one can see a sharper drop for the larger system 
sizes. The size dependence at U = 4, V = 2.1 is shown 
in an inset. Here the convergence to 0, i.e., the presence 
of a spin gap, is apparent. If the spin gap is small, as it 
is close to the phase boundary, the convergence to will 
obviously occur only for very large systems. 

Results such as those shown in Fig. [3 were used to de- 
termine the phase boundaries in Fig.^ As already noted, 
the BOW aspect of the phase diagram differs from pre- 
vious proposals [IS [H El El in that the BOW-CDW 
transition can be either continuous or first-order, i.e., the 
change of order occurs on the BOW-CDW boundary. The 
existence of two special points, one where the transition- 
order changes and one at higher U where the BOW van- 
ishes, was also suggested by Jeckelmann |l4|. who, how- 
ever, insists that the BOW does not exist for small U 
where the transition to the CDW state is continuous (i.e., 
his phase diagram has no continuous BOW-CDW tran- 
sition). We have carried out calculations for U down to 
1, and, as shown in Fig. ^ we still find a BOW phase 
there. Most likely, in view also of weak-coupling argu- 
ments 0, the BOW extends down to U = 0+ We find 
no BOW for U > 9. In the strong-coupling expansion, 
using the couplings J and J' derived by van Dongen [8j, 
the effective spin model is gapped, i.e., J' > 0.241J pL7| . 
above the dashed curve in Fig. ^ The J — J' mapping is 
not applicable beyond the transition into the CDW state, 
which (the solid curve in Fig.QJ was previously calculated 
by comparing the 4th-order CDW and J — J' energies 
[8j. The 4th-order BOW region ends at U w 7, where 
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FIG. 3: QMC results for Luttinger charge exponent on the 
BOW-CDW boundary (solid circles). The inset shows the 
finite-size scaling for [7 = 4. 



the spin-gap curve crosses the CDW-transition curve. 
This is slightly lower than what we find based on QMC. 
The strong-coupling BOW extends down to smaller U, V, 
but clearly the 4th-order result cannot be expected to 
be quantitatively accurate in this region. Nevertheless, 
the spin-frustration mechanism consistently explains the 
presence of an SDW-BOW transition and an extended 
BOW phase. Spin-frustration was previously cited by 
Jeckelmann |l4j |. but, surprisingly, he used it in support 
of a BOW of vanishing extent. 

Next, we consider the nature of the BOW-CDW tran- 
sition. As discussed in Refs. [lcllliij. the continuous crit- 
ical point for small U is described by a Gaussian free 
(charge) boson theory, characterized by the parameter 
K p . At generic values of (repulsive) U, V, the leading 
"4fcp" umklapp process is present, and has scaling di- 
mension A4k F = 2K p and is hence relevant (A^ < 2) 
for K p < 1. At the BOW-CDW transition, this oper- 
ator vanishes, leading to a vanishing of the charge gap. 
For consistency, no other relevant operators should be 
present, which would otherwise require fine tuning to 
zero, making the Gaussian theory a multicritical point. 
The most dangerous candidate is the "8/ci?" umklapp pro- 
cess, with A 8 fe F = 4A 4j t F = 8K p , so a continuous Gaus- 
sian critical point is possible only for 1 > K p > 1/4. Even 
in this range, the Gaussian theory is an unusual critical 
point with non-universal behavior, e.g., the correlation 
length exponent v = 1/(2 — 2K p ). 

Extrapolated QMC results for K p on the BOW-CDW 
boundary are shown in Fig. |3| The finite-size corrections 
appear to be of the form 1/L a , with an U dependent 
exponent a. At U = 4, a » 1, as shown in the inset of 
Fig. 13 For larger U, a decreases rapidly and is difficult 
to determine accurately for U ^ 5. The extrapolated 
K p value at U = 5 in Fig. should be regarded as an 
upper bound. At U = 6, the extrapolated K p < 1/4, and 
hence we expect an eventual drop to 0. This is consistent 
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FIG. 4: Finite-size scaling of the BOW and CDW susceptibil- 
ities for U = 3 (left) and 4 (right). The symbols correspond to 
different system sizes as in Fig. |5] The dashed lines indicate 
the independently determined critical points. 



with clear signals of a first-order transition . Also at 
U = 5.5 there are signs of first-order behavior, e.g., in 



order parameter histograms such as those considered in 
Ref. [llj We believe that the change from a continuous 
to a first-order transition occurs between [7 = 5 and 5.5. 

What is the nature of the tricritical point at which 
the transition becomes first-order? The simplest scenario 
is that this is the last (marginally) stable point of the 
Gaussian fixed line, i.e. with K p = 1/4. This hypothesis 
predicts that the critical K p continuously approaches 1 /4 
as the tricritical point U = Ut is approached from below, 
as K p — 1 /4 ~ yj (U t — U)/Ut- We do not have sufficient 
data to verify this form, but it is consistent with a sharp 
drop to between U = 5 and 5.5 (Fig. required since 
at U = 5.5 the transition should be first-order. Hence, 
we favor this behavior over the a priori consistent (but 
less simple) possibility of a non-trivial "strong coupling" 
tricritical fixed point far from the Gaussian line. 

To further demonstrate the Luther-Emery state on the 
continuous BOW-CDW curve, we study the finite-size 
scaling of the CDW and BOW susceptibilities, xcdw("") 
and xbow(tt) (with their standard Kubo- integral defini- 
tions Both the charge and bond correlations should 
decay as (— iyr~ K e 0, implying that the susceptibili- 
ties scale with system size as L 2 *>. Thus, x( 7T )L Kp ~ 2 
curves for different L should intersect at the critical 
BOW-CDW point. Fig. g| shows results for U = 3 and 
4, using the K p values determined above. For U — 4 the 
expected scaling can be observed even for small systems. 
For U — 3 the corrections are larger, and the asymptotic 
scaling sets in only for L > 128. This is clearly due to the 
smaller spin gap at U — 3, which implies a longer length- 
scale below which remaining spin correlations affect the 
charge and bond fluctuations. 

In summary, we have determined the ground state 



phase diagram of the extended Hubbard model at half- 
filling. The dimerized BOW phase can be explained 
by spin-frustration. The BOW-CDW transition changes 
from continuous to first-order between [7 = 5 and 5.5. 
On the critical (U, V) curve the system is a Luther- Emery 
liquid, with a charge exponent K p decreasing from 1 as U 
is increased from 0. We have argued that the minimum 
K p = 1/4 and that the BOW-CDW transition becomes 
first-order when this value is reached. 
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